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Abstract 

We investigate the position of the radio eore in a blazar by multi-epoch astrometric 
observations at 43 GHz. Using the VLBI Exploration of Radio Astrometry (VERA), we 
have eonducted four adjacent observations in Eebruary 2011 and another four in October 
2011, and succeeded in measuring the position of the radio core in the TeV blazar Mrk 501 
relative to a distant eompact quasar NRAO 512. During our observations, we find that (1) 
there is no positional ehange within ~ 0.2 mas or ~ 2.0 pc de-projected with ±la error for 
the weighted-mean phase-referenced positions of Mrk 501 eore relative to NRAO 512 over 
four adjaeent days, and (2) there is an indieation of position ehange for 3C 345 eore relative 
to NRAO 512. By applying our results to the standard internal shoek model for blazars, we 
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constrain the bulk Lorenz factors of the ejecta. 

Key words: astrometry — phase-referencing — galaxies: active — galaxies: blazars: 
individual (Mrk 501 = DA 426 = J1653-I-397; 3C 345; NRAO 512) — radio continuum: 
galaxies — techniques: interferometric 


1. Introduction 

Locating radio-emitting regions in relativistic jets in active galactic nuclei (AGNs) has been 
one of the most intriguing issues for exploring the ultimate mechanism of jet formation (e.g. Marscher 
et ah, 2008; Abdo et ah, 2010). One of the leading scenarios suggests that radio cores of blazars 
correspond to stationary standing shocks located at several parsec scales downstream from the central 
black holes, based on the delay timescales of the radio core brightening after the precedent y-ray 
flares (Marscher et ah, 2008; Jorstad et ah, 2010; Agudo et ah, 201 la,b). On the other hand, phase- 
referencing Very Long Baseline Array (VLBA) observations of the radio galaxy M87 revealed that 
the radio core at 43 GHz is located only ~ 0.01 pc away from the upstream end of the conical jet 
(Hada et ah, 2011). Nowadays, it is important to explore the origin of such a discrepancy of radio 
core locations and their stationarity between radio galaxies and blazars. 

To explore the locations and their stationarity of radio cores in blazars directly, multi-epoch 
Very Long Baseline Interferometry (VLBI) astrometric observations are the most powerful and reli¬ 
able method. By means of astrometric observations, we can measure the positions of target sources 
relative to phase-calibrators close to the targets on the sky. Over the past few decades, astrometric 
monitoring experiments based on the International VLBI Service for Geodesy and Astrometry (IVS) 
also have been performed towards the radio cores of more than a hundred of blazars (e.g., Ma et ah, 
1998; Fey et ah, 2004; Titov et ah, 2011). However, most of these astrometric experiments have been 
limited at 2.3 & 8 GHz'. Hence, little is known about the relative location of brightness peak and the 
stationarity seen at higher frequencies. 

We should carefully select an appropriate blazar and an observing frequency. It is well known 
that radio cores of blazars usually tend to be optically thick against synchrotron self-absorption (SSA) 
at lower frequencies (e.g., Kellermann & Pauliny-Toth, 1981; Marscher & Gear, 1985). To avoid 
the SSA effect as much as we can, we select a TeV blazar, Mrk 501 {z = 0.034). Mrk 501 is one 
of the best blazars for this study because the radio core spectrum between 1.4 GHz and 86 GHz 
measured by VLBI clearly shows convex shape peaking at ~ 8 GHz (figure 8 in Giroletti et ah, 
2008). Therefore, we can probe the locations of radio emitting regions reducing the SSA effect by 
selecting the observing frequency higher than the peak frequency. We select the observing frequency 
at 43 GHz because of the transparency against SSA. The proximity and brightness of Mrk 501 enable 
us to perform the precise astrometric observations. Another advantage of this source is that there are 

' see however Lanyi et al. (2010) 


2 



two bright phase-calibrators within an angular separation of a few degrees, such that it makes easier 
to perform phase-referencing observations. 

The plan of this paper is as follows. In §2, we describe our astrometric procedure for measuring 
the core position of Mrk 501. We explain the observation in §3, and data reduction in §4. The results, 
and the discussions are described in §5 and §6. The conclusion and future prospects are given in 
§7. The error estimation of the phase-referenced positions is written in appendix 1. Throughout this 
paper, we adopt the following cosmological parameters; Hq = 71 kms“^Mpc“\ ilu = 0.27, and 
IIa = 0.73 (Komatsu et ah, 2009), or 1 mas = 0.662 pc for Mrk 501. 

2. Astrometric procedure in measuring Mrk 501’s core position 

In this section, we describe our method to measure the core position of the target source 
Mrk 501. We chose the distant point source NRAO 512 (z = 1.66, 1 mas = 8.56 pc) as a reference 
position, since the absence of the jet component minimize the positional change of its location. The 
angular separation between Mrk 501 and NRAO 512 is 2?57. Given the angular separation limit of 
dual-beam system for VERA is 2?20, we inserted a bright calibrator 3C 345 (z = 0.593, 1 mas = 
6.64 pc), separated by 2.°09 from Mrk 501. The source configuration on the celestial sphere is shown 
in figure 1. Since the flux density of 3C 345 (~ 2 Jy) is much higher than that of Mrk 501 
300 mJy) and NRAO 512 500 mJy), we chose 3C 345 as a phase calibrator and applied its fringe 

phase solutions and structure phase solutions to the other two sources. The two sources, 3C 345 and 
NRAO 512, are a well-known phase-referencing pair (Shapiro et ah, 1979; Bartel et ah, 1986; Doi 
et ah, 2006; Jung et ah, 2011). Below, we explain our steps to derive the core position of Mrk 501 
relative to NRAO 512 for each epoch. 

(1) We measured the core position offset from the phase-tracking center on the phase- 
referenced images of Mrk 501 referenced to the 3C 345 core. 

(2) We measured the peak position offset from the phase-tracking center on the phase- 
referenced images of NRAO 512 relative to 3C 345 core. 

(3) We subtracted the position offset of (2) from that of (1) to derive the core position of 
Mrk 501 relative to NRAO 512. 

Through this procedure, the position errors arise from 3C 345 were completely cancelled out 
(for details, see §5.3 and A. 1.3). 

3. Observations 

The VERA observations were carried out for two sets of four adjacent days, on Eebruary 15, 
16, 17, 18, and October 20, 21, 23, 24, 2011 (details shown in table 1). Two sources within 2?2 
were observed simultaneously by the VERA’s dual-beam system (Honma et ah, 2003). In the total 
bandwidth of 256 MHz (16 X 16 MHz), one of the sixteen 16-MHz intermediate frequency (IE) 
channels was assigned to the bright calibrator, 3C 345 in A-beam. The other 15 IE channels were 
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allocated to the target Mrk 501, and NRAO 512 in B-beam (see also figure 1). Left hand circular 
polarization (LHCP) signals were received in the Q-band, 43115-43131 MHz in A-beam and 42987- 
43227 MHz in B-beam. The signals were digitized with 2-bit quantization by using the VERA digital 
filter (Iguchi et ah, 2005), and recorded with a rate of 1024 Mbps. 

The total observation time is 8 hours per each epoch, in which the total on source time is 
~ 8 hours for 3C 345 in A-beam, and ~6 hours for Mrk 501, ~2 hours for NRAO 512 in B-beam. 
The real-time instrumental phase difference data between the two beams were measured with artificial 
noise sources (Kawaguchi et ah, 2000; Honma et ah, 2008b). Data correlation was performed with 
the Mitaka FX correlator (Chikada et ah, 1991). 

4. Data reduction 

Initial calibration on correlated visibilities was performed with the Astronomical Image 
Processing System (AIPS) software package developed by the National Radio Astronomy 
Observatory (NRAO). At first, the visibilities were normalized by the auto-correlation with the 
AIPS task ACCOR and BPASS. A priori amplitude calibration was performed with the AIPS 
task APCAL on the basis of the system equivalent flux density (SEED) derived from the opacity- 
corrected system temperature and the antenna gain information for each antenna. The opacity- 
corrected system temperature was measured with chopper-wheel method. We adopted the antenna 
gain with a correction of its dependency on the separation angle between two beams, based on 
the VERA status report in 2009^. The accuracy of the amplitude calibration is ~10% (e.g., Nagai 
et ah, 2013; Petrov et ah, 2012). Secondly, we applied the recalculated delay-tracking solutions 
to the correlated data (Honma et ah, 2007). We adopted the delay-tracking center as (q:j 2 ooo> 
52ooo)=(16^53“52f216685, +39°45'36"60894) for Mrk 501, (aj2ooo, 52ooo)=(16M2“58."809965, 
+39°48'36"99399) for 3C 345, and (aj2ooo, ()2ooo)=(16MO“39f632772, +39°46'46"02849) for 
NRAO 512, taken from the International Celestial Reference Frame 2 (ICRF2, Fey et al. 2010). 
We also applied instrumental delay correction table to the data sets on B-beam in order to calibrate 
dual-beam delay difference (Honma et ah, 2008b). After that, we performed fringe fitting on the 
calibrator 3C 345 in the AIPS task FRING to solve delays, rates, and phases by assuming a point 
source model. The residual phase and gain of 3C 345 were solved by the self-calibration imaging 
(described in the next paragraph). With the AIPS task CAEIB, the phase and gain solutions of self¬ 
calibration were obtained using the source structure model. Thirdly, we transferred the derived phase 
solutions of 3C 345 to NRAO 512 and Mrk 501. After applying the amplitude solutions of Mrk 501 
and NRAO 512 obtained by the self-calibration imaging, we obtain the phase-referenced images of 
these two sources with the CEEAN procedure in Difmap software package (Shepherd, 1997). 

We constructed self-calibrated images of all sources. To perform amplitude self-calibration 
after phase self-calibration converged, we need all the four stations. However, since part of our 

^ http://veraserver.mtk.nao.ac.jp/restricted/CFP2009/status09.pdf 
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observations was lacking in one or two stations due to bad weather conditions (see table 1), we 
combined the visibilities over four adjacent observations by AIPS task DBCON to increase the tin- 
coverage for obtaining precise images. Imaging was performed in Difmap. After careful flagging of 
bad visibilities and averaging over ~5 seconds in the time domain, we constructed the initial structure 
models by model fitting and self-calibration iteratively. The self-calibrated image parameters are 
summarized in table 2. 

The peak position offsets from the phase-tracking center on the phase-referenced images of 
Mrk 501 and NRAO 512 were measured with the AIPS task JMFIT. The peak positions measured 
by model-fitting to the phase-referenced visibilities are consistent with those by JMFIT (typically 
within ~ 10 /ras, less than 30 fxas) and did not change the results significantly. We derived their 
phase-referenced core positions by subtracting from the JMFIT peak positions of the phase-referenced 
images the peak position differences between JMFIT to the self-calibrated images and model-fitting 
to the self-calibrated visibilities (~ 1 /ras). Then we subtracted the core position of NRAO 512 
referencing to 3C 345 from that of Mrk 501 referencing to 3C 345. Thus, we finally obtained the core 
position of Mrk 501 relative to NRAO 512. 

5. Results 

5.1. Phase-referenced image qualities 

Figure 2 shows the phase-referenced images of Mrk 501 relative to 3C 345. The phase- 
referenced images appear a point-like source due to the low SNR of the data. We also find NRAO 512 
relative to 3C 345 in the phase-referenced images as a single-component source. The image qualities 
of all the phase-referenced images are summarized in column (l)-(5) of table 3. We detect both of 
Mrk 501 and NRAO 512 with an image SNR of higher than nine for all epochs except for the third 
one, in which two stations suffered from bad weather condition over more than half of the observation 
time. By comparing the peak fluxes of the phase-referenced images in column (3) of table 3 to those 
of self-calibrated images in column (2) of table 2, the averaged flux loss ratio is evaluated as typically 
~ 20%, up to 60% under bad weather conditions. 

In appendix I, the procedure to estimate the position errors is described, and the error values 
are summarized in table 4. The position errors are dominated by the tropospheric zenith delay error, 
by assuming typically ~2 cm error for VERA (Honma et ah, 2008a). Since theoretically estimated 
thermal noise error (a few /ras, where the image sensitivity is estimated as 1 mJy beam“^ by as¬ 
suming the system noise temperature ~ 300 K) does not contribute to the random error (~ 20 /ras) so 
much, some other random errors and calibration errors still affect the quality of the phase-referenced 
images (e.g., Doi et ah, 2006). Here we address them as possible as we can. Firstly, there would 
be random fluctuation of tropospheric error and geometrical error. Assuming the typical fluctuation 
of the propagation delay as ~ 1 cm day“^ (Treuhaft & Lanyi, 1987), the position error would be 
~ 60 /ias for Mrk 501-3C 345 pair, and ~ 14 /ras for NRAO 512-3C 345 pair. Secondly, the errors of 
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self-calibration for 3C 345 are also related to the phase-referenced image quality; an antenna-based 
phase solution with a SNR of ~ 8 corresponds to an accuracy of 7° (~ 12 /ras) theoretically. Thirdly, 
there would be some errors when we identify the core positions of Mrk 501 and NRAO 512 in the 
phase-referenced images with Gaussian fitting, because they are not point sources. We evaluate the 
core identification error in the phase-referenced images as ~ 5 /ras, by summing the core identifi¬ 
cation errors (~a few /ras in A. 1.3) and the Gaussian fit errors (~a few /ras) in the self-calibrated 
images. However, for most epochs, the image SNR of NRAO 512 is comparable to that of Mrk 501, 
although the angular separation between NRAO 512 and 3C 345 is one-fourth of that of Mrk 501 and 
3C 345. This means that the image SNR would not be mainly limited by the error terms depending 
on the angular separations. We would conservatively overestimate the systematic error (especially for 
the tropospheric zenith delay error). 

5.2. Mrk 501 core position relative to NRAO 512 

Figure 3 shows the measured core positions of Mrk 501 relative to NRAO 512 for all the ob¬ 
servations as red points and blue points, which are listed in table 5. The core position error of Mrk 501 
relative to NRAO 512 for each epoch in table 5 is estimated by adding the root sum square of the ran¬ 
dom error (cirandom) to all the other systematic position errors of Mrk 501 and NRAO 512 relative 
to 3C 345 in table 4, except for the core identification error of 3C 345 (see §5.3 and appendix 1). 
Typically, the core position error of Mrk 501 relative to NRAO 512 for each epoch is ~ 0.20 mas in 
right ascension (RA) and ~ 0.21 mas in declination (Dec), which is around one-third of the major 
axis of the beam size. 

The core positions of Mrk 501 relative to NRAO 512 are distributed clustering during four 
adjacent days in February and October 2011 in the self-calibrated image. During the four adjacent 
days, the random error, the tropospheric error, and the ionospheric error (arandom, (^trop, and aion in 
table 4) vary randomly. All the other errors, which are the earth orientation parameter error, the 
antenna position error, a priori source coordinate error, and the core identification error (Uearth, o-ant, 
O'coordj and did in table 4) are systematic errors. The maximum differences of the core positions within 
the four days are 0.11±0.18 mas and 0.08 ±0.18 mas in February and October 2011, respectively. 
For the errors of the maximum position differences, the common systematic errors are canceled out. 
We perform a chi-squared test of the phase-referenced core positions over the four adjacent days, and 
confirm that the core positions of Mrk 501 relative to NRAO 512 coincide in both RA and Dec with 
the significance probability of > 98%. The weighted mean positions of Mrk 50rs core relative to its 
phase-tracking center over the four adjacent days are (x,y) = (—0.03 ± 0.10, —0.05 ± 0.11) mas in 
February 2011, and (x,y) = (—0.09 ± 0.10, 0.01 ± 0.10) mas in October 2011, shown as black points 
in figure 3. Both of random errors and systematic errors of each epoch are used as the weight for 
calculating the weighted mean positions and their conservative errors. The spatial distribution ranges 
of the core positions over the four adjacent days are estimated to be 0 ~ 0.22 mas and 0 ~ 0.20 mas 
in February and October 2011, respectively, with the statistical significance of ±1 a (68% confidence 
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level) for the weighted mean position errors. The eore positions of Mrk 501 relative to NRAO 512 
between February and Oetober 2011 (red erosses and blue erosses in figure 3) eould be systematieally 
different, but the differenee between the eenters of the weighted-mean positions is ~ 0.06 mas in 
both RA and Deo direotion, whieh is less than the spatial distribution ranges of the four adjaeent days 
0.20 mas). Therefore, we find there is no positional ehange of the radio eore in Mrk 501 relative 
to NRAO 512 within ~ 0.20 mas. The disoussion is shown in §6.1. 

5.3. Core positions of Mrk 501 and NRAO 512 relative to 3C 345 core 

Figures 4 (a) and (b) show the eore positions of Mrk 501 and NRAO 512 relative to 3C 345 
eore, respeetively. The phase-refereneed image qualities, the eore positions and their position errors 
are summarized in table 3. In the previous subseetion, we obtain the eore position of Mrk 501 relative 
to NRAO 512 for eaeh epoeh by subtraeting the eore position in figure 4 (b) from that in figure 4 (a). 
Contrary to the derived position offsets from the phase traeking eenter of Mrk 501 eore relative to 
NRAO 512 (within ~0.09 mas), the position offsets from the phase tracking center of Mrk 501 core 
relative to 3C 345 core and NRAO 512 relative to 3C 345 core are ~ 0.4 mas in RA and ~ 0.2 mas in 
Dec, which mainly caused by the position shifts of 3C 345 core. 

5.3.1. Mrk 501 core position relative to 3C 345 core 

The position errors of Mrk 501 core relative to 3C 345 core are ~ 0.12 mas in RA and ~ 
0.13 mas in Dec, summarized in column (15) and (16) of table 4. From figure 4 (a), the core positions 
of Mrk 501 relative to 3C 345 distributes clustering over the four adjacent days. The weighted mean 
positions are {x,y) = (—0.34 ± 0.06,—0.19 ± 0.06) mas in February 2011, and {x,y) = (—0.45 ± 
0.06,—0.20±0.06) mas in October 2011. The weighted mean position in February 2011 is ~0.11 mas 
east of that in October 2011, which is less than the distribution range (~ 0.13 mas with ±la error). 
Thus, there is no positional change of Mrk 501 core relative to 3C 345 core within ~ 0.13 mas. 
Note that the phase-referenced positions may include the possible positional change of Mrk 501 
« 0.20 mas) and 3C 345, because both of them are not point sources. Therefore, we do not discuss 

this result in further detail. 

5.3.2. NRAO 512 position relative to 3C 345 core 

The position errors of NRAO 512 relative to 3C 345 core are ~ 0.04 mas both in RA and 
Dec (table 4). The estimated position errors are about one-fourth of those for Mrk 501 core relative 
to 3C 345 core mainly due to the tropospheric error based on one-fourth smaller angular separation. 
In figure 4 (b), we find the core positions of NRAO 512 relative to 3C 345 in February 2011 (red 
points) and October 2011 (blue points) are clearly different. The weighted mean positions are (x,?/) = 
(-0.31 ± 0.02, -0.14 ± 0.02) mas in February 2011, and {x,y) = (-0.36 ± 0.02, -0.21 ± 0.02) mas 
in October 2011. The weighted mean position in October 2011 is 0.05 ± 0.03 mas west and 0.07 ± 
0.03 mas south of that in February 2011. This result indicates peak position change of 3C 345 by 
assuming the peak position of NRAO 512 is stationary. The discussion is shown in §6.2. 
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6. Discussions 


6.1. Application to internal shock model 


Here we attempt to constrain bulk Lorentz factors of the jet based the result shown above. We 
assume that the location of the unresolved radio emitting regions in Mrk SOl’s core at 43 GHz are 
identical to that of the soft X-ray emitting regions since the radio core is optically thin against SSA. 
Therefore, the internal shock model can be applicable to our target source Mrk 501.^ The internal 
shock model can well explain the various observational properties of soft X- and gamma-ray light- 
curves in blazars and thus is regarded as one of the leading and standard models for blazar emissions 
(e.g., Spada et ah, 2001; Tanihata et ah, 2003; Guetta et ah, 2004; Kino et ah, 2004; Mimica et ah, 
2004; Bottcher & Dermer, 2010; Joshi & Bottcher, 2011). Internal shocks will occur when a faster 
discrete ejecta catches up with slower one. Electrons are accelerated by the shocks and the non- 
thermal electromagnetic waves are emitted by these relativistic electrons (e.g., Mimica & Aloy, 2010, 
2012). Then, denoting the faster ejecta’s Lorentz factor (Tf), slower ejecta’s Lorentz factor (Tg), and 
the initial separation of the ejecta (/is), the de-projected distance between the internal-shock position 
and the central engine (T*is) is given by: 


r?r2 


D 


IS 


fs 


'rf-r? 


/is = 2- 


(ri/r. 


-r?/is 


(rf/r,)2-i ^ 

where Tf > Tg 3> 1. In general, such collisions happen repeatedly and multiple internal-shocked 
regions are generated in the jet flow. The de-projected distributed scale of the shocked regions is, by 
definition, given by 


ADis = D 


IS,max 


D 


IS,min? 


( 1 ) 


where i/is,max is the largest distance between the location of the internal-shocked region and that of 
the central engine, and i/is,min is the closest one. 

The term ADi^ for Mrk 501 is directly constrained by our VERA astrometric observations. 
In §5.2, we reveals that the weighted mean core positions of Mrk 501 over the four adjacent days 
in Lebruary and October 2011 were spatially distributed within 0 ~ 0.22 mas and 0 ~ 0.20 mas, 
respectively, by assuming the radio core of NRAO 512 is stationary. The de-projected distributed 
scale ADis can be expressed as follows by using the lower limit of the jet viewing angle of 0j>4° 
(Giroletti et ah, 2004): 


AL)is < 2.2 X 10^ /?s 


0 


0.20 mas 



( 2 ) 


^ Here we do not apply our result to the standing shock model. Recently the original model of the conical standing 
shock proposed by Marscher et al. (2008) has been improved considering the multi-zone in the standing shock region 
(so called ’’turbulent extreme multi-zone” (TEMZ) model) by Marscher (2014), which also can explain the variability 
of the soft X-ray light curve of the blazar. For now, our result does not rule out the model because the position of the 
standing shock is stationary. 
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or ADis <1.9 pc, where 1 mas corresponds to 7.7 x 10^ Schwarzschild radii (i?s), or 1 i?s = 8.6 x 10“^ 
pc. Here we set the lower limit of the central black hole mass of Mrk 501 as Mbh = 0.9 x 10® Mq in 
the case of the single black hole (Barth et ah, 2002). 

The right-hand-side of Eq. (1) has been constrained by the previous work. We set the sepa¬ 
ration between the ejecta satisfies Ijs > IRs, because IRs is the minimum dimension of the central 
engine (e.g., Spada et ah, 2001). Tanihata et al. (2003) suggested that the Lorenz factor ratio between 
the faster ejecta and the slower ejecta is Ef/Es < 1.01 to reproduce the observed X-ray light curves 
in the flare state. In quiescent state, Ef/Es would be less than that in the flare state because larger 
Ef/Es produce stronger internal shocks closer to the central engine and generate larger energy than 
smaller Ef/Es do when we assume constant dynamical efficiency and constant Jfs. Therefore, we can 
assume that (Ef/Es)Dis=Dis (^i/^s)Dis=Dismin ^ 1-01 because during our observation Mrk 501 
was in relatively quiescent state (just before the X- and y-ray flare; Bartoli et al. 2012). When we 
assume Ef/Eg is constant, His,max and His,min only depend on the maximum value of Eg (Eg,max) and 
minimum value of Eg (Eg,min), respectively. Now, AHis can be expressed as a function of only Eg. 
Regarding the minimum value of Eg, we adopt Eg min > 8 in Mrk 501 (e.g.. Kino et al., 2002) based 
on the minimum among the references of the one-zone SED model fitting. Using Eqs. (1) and (2), we 
finally constrain a maximum of Eg max/Eg,min by applying the maximum or minimum values of above 
assumptions as follows: 


s.max ^ o 1 

-< 2.lx 


s,min 


AH 


IS 


2.2 X lO^ H, 


-i 




77l\ 


V51/ VIH 


-1 


-2 


-,1/2 


+ 1 


(3) 


where A = (Ef/Es)^/[(Ef/Eg)^ — 1] and A > 51. Thus, we find that the maximum-to-minimumratio of 
Eg during our four adjacent observations are less than 2.1, i.e., 8 < Eg < 17, including the uncertainties 
of the relevant quantities. The derived Eg is roughly comparable to the one estimated by broadband 
spectra during its quiescent state (e.g., Abdo et al., 2011). Eurther observations during its flare state 
are encouraged to obtain larger AHis. 

Using the above Eg, we estimate His, which corresponds to the distance between the location 
of the central engine and that of the internal shock, where corresponds to the location of the radio core. 
The estimated His are His,min ~ 0.6 pc , and His,max ~ 2.5 pc (, and they are comparable 
to the distance suggested in other blazars (Marscher et al., 2008; Jorstad et ah, 2010; Agudo et al., 
2011a,b). 


6.2. Indication of peak position shift between NRAO 512 and 3C 345 

Despite the non-detection of the core position change of the target source Mrk 501, we find 
a shift of the relative peak positions between the two calibrators, NRAO 512 and 3C 345 (§5.3.2). 
The relative peak positions of the two quasars are clustering over four adjacent days in Eebruary and 
October 2011, respectively, and the positions of the clusters are different (see red points and blue 
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points in figure 4-b). From §5.3.2, the weighted mean position of NRAO 512 relative to 3C 345 core 
in October 2011 is 0.05 ± 0.03 mas west and 0.07 ± 0.03 mas south of that in February 2011. This 
may indicate a shift of 3C 345 core of 0.05 ± 0.03 mas east and 0.07 ± 0.03 mas north from February 
to October 2011, when we assume NRAO 512 is a stationary point source. The proper motion of 
3C 345 core is estimated to be 0.12 ± 0.04 mas yr“^, equivalent to apparent velocity 4.2 ± 1.4c at a 
position angle of 36^27 deg (measured from north to east). Bartel et al. (1986) performed astrometry 
between this phase-referencing pair mainly at 8 GHz, and find the proper motions of the core are 
not significantly different from zero with weighted least-squares fit to the nine epoch data spanned 
over ~ 11 years. However, focusing on their first three epoch data within one year, the positional 
change of 3C 345 core was observed toward similar direction. As for the 3C 345 jet ~ 0.60 mas 
away from the core, we find the proper motion of the absolute position ~ 0.40 mas yr“^ (aparent 
velocity ~ 14c), which is almost consistent with the result in Bartel et al. (1986), toward at a position 
angle of about —98° (or 72°). Although there is a small difference in the position angle of the proper 
motion between the jet and the core, a shift of the 3C 345 core could be related to the formation of 
the new jet component, because the large variation of the jet ejection angle has been observed with an 
amplitude of ~ 20° (e.g., Klare et al., 2005; Schinzel et al., 2010). This encourages further astrometric 
observations. 

7. Conclusion 

In order to investigate the stationarity of the radio-core position in blazars, we have conducted 
multi-epoch astrometric VERA observations of TeV blazar Mrk 501 at 43 GHz for the first time. 
Below, we summarize the main conclusions in the present work. 

1. In our observation, we achieve an accuracy of ~0.20 mas for the radio core position of Mrk 501 
relative to that of NRAO 512 for each epoch. The relative radio-core positions of Mrk 501 over 
four adjacent days in February 2011 are distributed within ~ 0.22 mas spatial scale, and those 
in October 2011 are within ~ 0.20 mas, respectively, with Tier error (68% confidence level) 
of the weighted-mean position. Comparing the weighted-mean position of the four adjacent 
observations in February 2011 to that obtained in October 2011, the difference between these 
two positions is ~ 0.07 mas. Summing up, we find that the radio-core position of Mrk 501 
during our observations does not show significant positional change within about 0.20 mas, by 
assuming the radio-core position of NRAO 512 is stationary. 

2. By assuming the standard internal shock model, we further constrain the bulk Lorentz factors of 
the ejecta based on the observational results, i.e., the core position stationarity within 0.20 mas 
(1.9 pc de-projected). Then we find that the maximum-to-minimum ratio of the slower ejecta’s 
Lorenz factor can be constrained to be rs,max/rs,min < 2.1. The distance between the location of 
the radio core and that of the central black hole is estimated as a few pc, which is similar to the 
one indicated in other blazars. 

3. An indication of core position change between February and October 2011 was found for the 
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phase-referencing pair of 3C 345-NRAO 512. This detection is quite encouraging for the future 
subsequent observations of blazar core astrometry initiated in the present work. 

Further continued astrometric observations targeting large X-ray flares of TeV blazar Mrk 501 are 
necessary to explore the locations and their stationarity of the radio core. Any simultaneous gamma or 
radio lightcurves are also interesting, as well as potential changes of the optical polarization position 
angles (e.g., Marscher et ah, 2008). 

The VERA is operated by Mizusawa VLBI Observatory, a branch of National Astronomical 
Observatory of Japan. We thank K. M. Shibata, T. like, and all the staff who helped operations of 
the VERA observations presented in this paper. We are grateful to T. Oyama, T. Hirota, M. Kim, N. 
Matsumoto, M. Sato for discussing how to reduce the VERA data. S.K. thanks the internal referee at 
the MPIfR, E. Ros, for his constructive comments. We are also grateful to S. Mineshige, K. Kohno, 
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24340042 (A.D.) and 2450240 (M.K.) from the Japan Society for the Promotion of Science (JSPS). 
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Mrk501 3C345 NRA0 512 

(B-beam) (A-beam) (B-beam) 



Fig. 1. Source configurations, dual-beam allocations, and self-calibrated images of our VERA astrometric observa¬ 
tions at 43 GHz. The three sources quasi-perfectly align within ~3' in east-west direction on the sky. The calibrator 
3C 345 was observed in A-beam, and at the same time, Mrk 501 and NRAO 512 were observed in B-beam in 
every several minutes. Shown are the self-calibrated images of the observed sources stacked over February 15, 
16, 17, 18, 2011. The contours are plotted at 3cr noise level x (-1, 1, 2, 4,). The size of the restored beam 
(bottom left) is 0.613 mas x 0.366 mas at a position angle at —53°. The Itr noise level is 3.04 mJy beam“^ 
for Mrk 501, and 5.15 mJy beam“^ for NRAO 512. The Icr noise level is 30.0 mJy beam“^ for 3C 345, 
which is about ten times higher than that for Mrk 501 because the image quality of the self-calibrated image 
is mainly limited by the image dynamic range (the peak flux of 3C 345 is about ten times higher than that of 
Mrk 501, see table 2). Dotted arrow shows the innermost jet direction of 3C 345 at a position angle of —103°. 
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Table 1. Summary of VERA 43 GHz astrometric observation 


Epoch 

time range (UT) 

Code 

bad weather 

2011/02/15 

18:18-02:18 

R11046B 

OGA: 0:30-1:20, 2:15-02:18, IRK: 0:50-1:40 

2011/02/16 

18:14-02:14 

R11047B 

IRK, OGA: 22:00-22:40 

2011/02/17 

18:10-02:10 

R11048B 

IRK, ISH: 18:00-23:00 

2011/02/18 

18:06-02:06 

R11049B 

IRK, OGA: 23:40-23:58,00:45-1:07 

2011/10/20 

01:50-09:50 

R11293A 

OGA: 01:50-06:00, IRK: 06:00-09:50 

2011/10/21 

01:46-09:46 

R11294A 


2011/10/23 

01:38-09:38 

R11296A 

ISH: 01:38-06:30,05:30-07:00 

2011/10/24 

01:34-09:34 

R11297A 



Notes. —MIZ: Mizusawa, IRK: Iriki, OGA: Ogasawara, ISH: Ishigaki. MIZ couldn’t observe Mrk 501 for ~50 min in 
every observation due to its elevation limit. 
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Table 2. Self-calibrated image parameters for 3C 345, Mrk 501, and NRAO 512. 


Source 

Epoch 

Ip 

^rms 

-^p/ f^rms 



(mjy beam“^) 

(mJy beam“^) 



(1) 

(2) 

(3) 

(4) 

3C 345 

2011/02 

2210 

30.0 

74 


2011/10 

2040 

17.0 

120 

Mrk 501 

2011/02 

280 

3.04 

92 


2011/10 

239 

1.78 

134 

NRAO 512 

2011/02 

467 

5.15 

92 


2011/10 

541 

1.94 

278 


Notes.—(1) Observing epoch, (2) peak intensity, (3) rms of image noise, (4) image dynamic range. 
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Fig. 2. Phase-referenced images of Mrk 501 core relative to 3C 345 core. The restored beam size is 0.669 mas 
X 0.391 mas at a position angle at —4019 and the contour start from three times 28.0 mJybeam”^. 
The contours are plotted 3cr noise level x (-1, 1.41, 2, 2.82, 4). Note that the image SNR of the 
third epoch is much lower than that of all the others due to the bad weather condition (see table 3). 
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Table 3. Image qualities and core positions of phase-referenced images. 


Phase-referenced pair 

Epoch 

4 

^rms 

-^p/ ^rms 

Aa 

A(5 



(mjy beam“^) 

(mJy beam”^) 


(mas) 

(mas) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Mrk501-3C 345 

2011/02/15 

259 

29.1 

9 

-0.33 ±0.13 

-0.21 ±0.13 


2011/02/16 

262 

21.0 

12 

-0.37 ±0.12 

-0.15 ±0.13 


2011/02/17 

224 

44.9 

5 

-0.38 ±0.14 

-0.18 ±0.14 


2011/02/18 

275 

30.8 

9 

-0.30 ±0.13 

-0.23 ± 0.13 


2011/10/20 

262 

15.5 

17 

-0.43 ±0.12 

-0.21 ±0.12 


2011/10/21 

210 

24.4 

9 

-0.44 ±0.13 

-0.19 ±0.13 


2011/10/23 

254 

18.5 

14 

-0.50 ±0.12 

-0.23±0.13 


2011/10/24 

270 

18.4 

15 

-0.44 ±0.12 

-0.15 ±0.12 

NRA0 512-3C 345 

2011/02/15 

388 

28.5 

14 

-0.33 ± 0.03 

-0.14 ±0.04 


2011/02/16 

408 

32.2 

13 

-0.31 ±0.03 

-0.15 ±0.04 


2011/02/17 

272 

61.6 

4 

-0.28 ± 0.08 

-0.12 ±0.07 


2011/02/18 

438 

27.5 

16 

-0.30 ± 0.03 

-0.14 ±0.04 


2011/10/20 

228 

17.4 

13 

-0.36 ± 0.04 

-0.18 ±0.04 


2011/10/21 

222 

25.2 

9 

-0.32 ± 0.04 

-0.23 ± 0.04 


2011/10/23 

424 

38.4 

11 

-0.38 ± 0.04 

-0.24 ± 0.04 


2011/10/24 

410 

32.5 

13 

-0.37 ± 0.04 

-0.18 ±0.04 


Notes.—(1) The first source is the target source and the second source is the reference source, (2) Observing epoch, (3) 
peak intensity of the target source, (4) image rms noise of the target source, (5) signal-to-noise ratio of the target source, 
(6) and (7) position offset from phase-tracking center of the target source. 
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Table 4. Estimated positional error budgets in the phase-referencing observation by VERA 43 GHz. The units are in /i.as. 


Phase-reference pair 

Epoch 

^ random 

fT’trop 

fT’ion 

^ earth 

<7ant 

^^coord 

CTid 


t^rss 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

Mrk501-3C 345 

2011/02/15 

27 

26 

122 

< 1 

3 

5 

1 

5 

2 

3 

7 

16 

125 

126 


2011/02/16 

16 

16 

122 

< 1 

3 

5 

1 

5 

2 

3 

7 

16 

123 

125 


2011/02/17 

56 

61 

122 

< 1 

3 

5 

1 

5 

2 

3 

7 

16 

135 

137 


2011/02/18 

27 

26 

122 

< 1 

3 

5 

1 

5 

2 

3 

7 

16 

125 

126 


2011/10/20 

13 

15 

122 

< 1 

3 

5 

1 

5 

2 

3 

9 

11 

123 

124 


2011/10/21 

23 

28 

122 

< 1 

3 

5 

1 

5 

2 

3 

9 

11 

125 

126 


2011/10/23 

16 

21 

122 

< 1 

3 

5 

1 

5 

2 

3 

9 

11 

123 

125 


2011/10/24 

16 

16 

122 

< 1 

3 

5 

1 

5 

2 

3 

9 

11 

123 

124 

NRA0 512-3C 345 

2011/02/15 

18 

16 

28 

< 1 

3 

5 

1 

5 

2 

3 

7 

16 

33 

38 


2011/02/16 

17 

17 

28 

< 1 

3 

5 

1 

5 

2 

3 

7 

16 

33 

37 


2011/02/17 

64 

70 

28 

< 1 

3 

5 

1 

5 

2 

3 

7 

16 

76 

72 


2011/02/18 

15 

14 

28 

< 1 

3 

5 

1 

5 

2 

3 

7 

16 

32 

36 


2011/10/20 

17 

18 

28 

< 1 

3 

5 

1 

5 

2 

3 

9 

11 

35 

36 


2011/10/21 

22 

29 

28 

< 1 

3 

5 

1 

5 

2 

3 

9 

11 

42 

38 


2011/10/23 

21 

24 

28 

< 1 

3 

5 

1 

5 

2 

3 

9 

11 

38 

38 


2011/10/24 

18 

20 

28 

< 1 

3 

5 

1 

5 

2 

3 

9 

11 

36 

36 


Notes.—(1) The first source is the target source and the second source is the reference source, (2) observing epoch, (3) and 
(4) the random errors estimated by the beamwidth over two times signal-to-noise ratio in RA and Dec, (5) the tropospheric 
residual errors, (6) the ionospheric residual errors, (7) and (8) the earth orientation parameter errors in RA and Dec, (9) 
and (10) the antenna position errors in RA and Dec, (11) and (12) a priori source coordinates errors in RA and Dec. The 
error contributions from the geometrical errors (7-12) are estimated based on the simulation presented in Pradel et al. 
(2006). (13) and (14) The core identification error of 3C 345 in RA and Dec, respectively. (15) and (16) Total errors in 
RA and Dec are estimated as the root-sum-square of each error. 
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Table 5. Core positions of Mrk 501 relative to NRAO 512. 


Epoch 

Act 

A() 


(mas) 

(mas) 

(1) 

(2) 

(3) 

2011/02/15 

-0.01 ± 0.20 

-0.10 ±0.20 

2011/02/16 

-0.06 ± 0.19 

0.00 ± 0.20 

2011/02/17 

-0.10 ±0.30 

-0.10 ±0.30 

2011/02/18 

0.01 ±0.19 

-0.10 ±0.20 

2011/10/20 

-0.07 ±0.19 

0.00 ± 0.20 

2011/10/21 

-0.10 ±0.20 

0.00 ± 0.20 

2011/10/23 

-0.12 ±0.19 

0.00 ± 0.20 

2011/10/24 

-0.07 ±0.19 

0.00 ± 0.20 


Notes.—(1) Observing epoch, (2) and (3) core position of Mrk 501 relative to NRAO 512 in RA and Dec, respectively. 
These values are derived by subtracting the core position of Mrk 501 relative to 3C 345 from tthat of NRAO 512 relative 
to 3C 345, summarized in (6) and (7) of table 3 for each epoch. Each position error is estimated as root-sum-square of 
corresponding error in table 3 (see §5.2). 
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Fig. 3. Self-calibrated image of Mrk 501 in figure 2 (a) is shown with the core positions of Mrk 501 relative to 
NRAO 512 for all epochs overlaid. The red points represent the core positions over four adjacent days in February 
2011 with an averaged positional accuracy of 0.21 mas in RA and 0.22 mas in Dec. The blue points show the core 
positions over four adjacent days in October 2011 with a positional accuracy of 0.19 mas in RA and 0.20 mas in 
Dec. Each core position is summarized in Table 5. The black points describe the weighted-mean positions over 
the four days in February 2011 at the relative (RA, Dec) = (—0.03 ± 0.10, —0.05 ± 0.11) mas, and (RA, Dec) = 
(—0.09 ±0.10, 0.01 ± 0.10) mas in October 2011. The origin of this hgure corresponds to the phase-tracking center 
of Mrk 501 and the brightness peak of the self-calibrated image. 



Fig. 4. (a) The core positions of Mrk 501 relative to 3C 345 core, and (b) those of NRAO 512 relative to 3C 345 
core. Red points represent four adjacent observations in February 2011, blue points represent in October 2011, and 
black points represent the weighted-mean positions over the four adjacent days (see §5.3). 
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Appendix 1. Error Estimation 

In this section, we estimate the position errors of the phase-referenced images referencing the 
error estimations adopted by previous VLBI astrometric studies (e.g., Hada et ah, 2011; Reid et ah, 
1999; Ros et ah, 1999). The error terms are the random errors and the uncertainties of ionospheric 
residuals, tropospheric residuals, core identification process, phase-connection process, instrumental 
origin, and geometrical parameters (earth orientation parameters, antenna positions and a priori source 
coordinates) (Thompson et ah, 2001). In the following subsections, details of major errors are de¬ 
scribed. Thanks to VERA’s dual-beam simultaneous phase-referencing system, the phase-connection 
errors can be ignored. The instrumental errors are canceled out by applying the dual-beam delay 
difference table measured with noise sources with ~ 0.1 mm accuracy (Honma et ah, 2008b). In 
table 4, we summarize the estimated error budget for each phase-referencing pair in our observations. 
We consider these errors contribute independently to each other, such that the total position errors for 
each epoch can be estimated as the root-sum-square of each error. 

A. 1.1. Random 

In column (3) and (4) of table 4, we summarize the random error of all phase-referenced 
images in RA and Dec, respectively. The position error can be expressed as 6*b/(2 x SNR) (Thompson 
et ah, 2001), where 9b is each interferometric beam size in RA or Dec (6*^ ~ 0.6 mas x 0.4 mas at a 
position angle of ~ 130°), and SNR is the signal-to-noise ratio of the phase-referenced images (see 
column 5 of table 3). The SNR of the phase-referenced images achieved is ~11 on average. 

A.1.2. Propagation delay 

Here we briefly review how to estimate the contribution of the propagation delay error to 
the phase-referenced position error. The effects of the propagation delay error at each station cause 
a decrease in SNR of the phase-referenced images and image shift. The position error 59 originated 
from the phase error A0 is expressed as 59 = A0 • (e-g-, Thompson et ah, 2001), where D\ ~ 

D ■ D is the baseline length, and A is the observing wavelength. AcJ) is estimated by 27rcArA“^, 
where c is the speed of light and At is the propagation delay error. At can be approximated by 
5to secZ, where 5to is the residual vertical delay and Z is the local source zenith angle. When 
the two phase-referencing pair separates in zenith angle by AZ, the position error 59 caused by the 
propagation delay error for a single antenna is estimated as a first-order Taylor expansion of r as 
follows (Reid et ah, 1999): 

59 ~ ggQ z tan Z AZ. (Al) 

The propagation delays are mainly caused by the ionospheric and tropospheric medium. In 
columns (5) and (6) of table 4, the tropospheric error and ionospheric error of each phase-referencing 
pair are summarized. The propagation delays caused by the tropospheric medium are independent of 
observing frequency, while those caused by the ionospheric medium are inversely proportional to the 
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square of the frequency. Therefore, we omitted the detail description of the ionospheric error, because 
it is estimated to be < 1 /ras at 43 GHz. Since the tropospheric zenith delays are typically within ~2 
cm accuracy for VERA (Honma et ah, 2008a), the position errors can be estimated as 66 ~ 122 /ras 
for AZ = 2.°09 (Mrk 501-3C 345 pair), and ~ 28 /ras for AZ = 0.°49 (NRAO 512-3C 345 pair) at 43 
GHz, Z ~ 50° and D = 2.3x 10^ m. 

A. 1.3. Core identification 

The core identification error indicates the possible slight difference between the fitted gaussian 
peak and the actual brightness peak, and would come from the blending of core and jet structure. 
To clarify the source structure, we produced the self-calibrated images of all the three sources (see 
figure 1 and §4). 

We identified the core of 3C 345 by model fitting to the calibrated visibility data, and found 
it has the bright jet structure. The core was well modeled by an elliptical Gaussian model, and the 
jet emission was modeled as one or two circular Gaussians. We defined the center of each Gaussian 
model as the component position. To evaluate the core identification error, we derived the differences 
between the center of the elliptical Gaussian and the position of the brightness-peak pixels in super¬ 
resolution maps measured by the AIPS task MAXFIT. The super-resolution maps were convolved by 
a circular Gaussian beam of full-width at half-maximum (FWHM) about a half of the minor axes of 
the synthesized beams. The maximum core identification error was estimated as 16 /ras, and all the 
position errors are summarized in column (13) and (14) of table 4. However, by subtracting the core 
position of NRAO 512 relative to 3C 345 from the core position of Mrk 501 relative to 3C 345, the 
core identification error of calibrator 3C 345 is completely cancelled out. 

For Mrk 501 and NRAO 512, we estimated the core identification error in the same man¬ 
ner. The position differences between visibility-based model-fitting and MAXFIT were negligible, 
typically less than a few /ras. 

References 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, Nature, 463, 919 
—. 2011, ApJ, 727, 129 

Agudo, L, Jorstad, S. G., Marscher, A. R, et al. 2011a, ApJF, 726, F13 

Agudo, L, Marscher, A. R, Jorstad, S. G., et al. 2011b, ApJF, 735, FIO 

Bartel, N., Herring, T. A., Ratner, M. L, Shapiro, I. L, & Corey, B. E. 1986, Nature, 319, 733 

Barth, A. J., Ho, F. C., & Sargent, W. F. W. 2002, ApJF, 566, F13 

Bartoli, B., Bemardini, R, Bi, X. J., et al. 2012, ApJ, 758, 2 

Bottcher, M., & Dermer, C. D. 2010, ApJ, 711, 445 

Chikada, Y., Kawaguchi, N., Inoue, M., et al. 1991, in Frontiers of VFBI, ed. H. Hirabayashi, 
M. Inoue, & H. Kobayashi, 79 
Doi, A., Fujisawa, K., Habe, A., et al. 2006, PASJ, 58, 777 


21 



Fey, A. L., Gordon, D., & Jacobs, C. S. 2010, in lERS/IVS Working Group, The Second Realization 
of the International Celestial Reference Frame by Very Long Baseline Interferometry, lERS 
Technical Note No. 35, http://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn35.html 
Eey, A. L., Ma, C., Arias, E. E, et al. 2004, AJ, 127, 3587 
Giroletti, M., Giovannini, G., Cotton, W. D., et al. 2008, A&A, 488, 905 
Giroletti, M., Giovannini, G., Taylor, G. B., & Ealomo, R. 2004, ApJ, 613, 752 
Guetta, D., Ghisellini, G., Lazzati, D., & Celotti, A. 2004, A&A, 421, 877 
Hada, K., Doi, A., Kino, M., et al. 2011, Nature, 477, 185 
Honma, M., Tamura, Y., & Reid, M. J. 2008a, PASJ, 60, 951 
Honma, M., Eujii, T., Hirota, T., et al. 2003, PASJ, 55, L57 
Honma, M., Bushimata, T., Choi, Y. K., et al. 2007, PASJ, 59, 889 
Honma, M., Kijima, M., Suda, H., et al. 2008b, PASJ, 60, 935 
Iguchi, S., Kkurayama, T., Kawaguchi, N., & Kawakami, K. 2005, PASJ, 57, 259 
Jorstad, S. G., Marscher, A. P, Larionov, V. M., et al. 2010, ApJ, 715, 362 
Joshi, M., & Bottcher, M. 2011, ApJ, 727, 21 
Jung, T., Sohn, B. W., Kobayashi, H., et al. 2011, PASJ, 63, 375 

Kawaguchi, N., Sasao, T., & Manabe, S. 2000, in Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, Vol. 4015, Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, ed. H. R. Butcher, 544-551 
Kellermann, K. L, & Pauliny-Toth, 1.1. K. 1981, ARA&A, 19, 373 
Kellermann, K. L, Vermeulen, R. C., Zensus, J. A., & Cohen, M. H. 1998, AJ, 115, 1295 
Kino, M., Mizuta, A., & Yamada, S. 2004, ApJ, 611, 1021 
Kino, M., Takahara, E, & Kusunose, M. 2002, ApJ, 564, 97 

Klare, J., Zensus, J. A., Lobanov, A. P, et al. 2005, in Astronomical Society of the Pacific Conference 
Series, Vol. 340, Euture Directions in High Resolution Astronomy, ed. J. Romney & M. Reid, 40 
Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 
Lanyi, G. E., Boboltz, D. A., Chariot, P, et al. 2010, AJ, 139, 1695 
Ma, C., Arias, E. E, Eubanks, T. M., et al. 1998, AJ, 116, 516 
Marscher, A. P. 2014, ApJ, 780, 87 

Marscher, A. P, Jorstad, S. G., D’Arcangelo, E D., et al. 2008, Nature, 452, 966 
Marscher, A. P, & Gear, W. K. 1985, ApJ, 298, 114 
Mimica, P, & Aloy, M. A. 2010, MNRAS, 401, 525 
—. 2012, MNRAS, 421, 2635 

Mimica, P, Aloy, M. A., Muller, E., & Brinkmann, W. 2004, A&A, 418, 947 
Nagai, H., Kino, M., Niinuma, K., et al. 2013, PASJ, 65, 24 
Petrov, L., Honma, M., & Shibata, S. M. 2012, AJ, 143, 35 
Pradel, N., Chariot, P, & Lestrade, J.-E. 2006, A&A, 452, 1099 

Reid, M. J., Readhead, A. C. S., Vermeulen, R. C., & Treuhaft, R. N. 1999, ApJ, 524, 816 


22 



Ros, E., Marcaide, J. M., Guirado, J. C., et al. 1999, A&A, 348, 381 

Schinzel, F. K., Lobanov, A. R, & Zensus, J. A. 2010, in Astronomical Society of the Pacific 
Conference Series, Vol. 427, Accretion and Ejection in AGN: a Global View, ed. L. Maraschi, 
G. Ghisellini, R. Della Ceca, & F. Tavecchio, 153 
Shapiro, 1. L, Wittels, J. J., Counselman, III, C. C., et al. 1979, AJ, 84, 1459 
Shepherd, M. C. 1997, in Astronomical Society of the Pacific Conference Series, Vol. 125, 
Astronomical Data Analysis Software and Systems VI, ed. G. Hunt & H. Payne, 77 
Spada, M., Ghisellini, G., Lazzati, D., & Celotti, A. 2001, MNRAS, 325, 1559 
Tanihata, C., Takahashi, T., Kataoka, J., & Madejski, G. M. 2003, ApJ, 584, 153 
Thompson, A. R., Moran, J. M., & Swenson, Jr., G. W. 2001, Interferometry and Synthesis in Radio 
Astronomy, 2nd Edition, 316, 467-506 
Titov, O., Lambert, S. B., & Gontier, A.-M. 2011, A&A, 529, A91 
Treuhaft, R. N., & Lanyi, G. E. 1987, Radio Science, 22, 251 


23 



